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We develop a bosonization approach to study the low tem- 
perature properties of one-dimensional gas of particles obey- 
ing fractional exclusion statistics (FES). It is shown that such 
ideal gas reproduces the low-energy excitations and asymp- 
totic exponents of a one-component Luttinger liquid (with no 
internal degrees of freedom). The bosonized effective theory 
at low energy (or temperature) is identified to a c = 1 confor- 
mal field theory (CFT) with compactified radius determined 
by the statistics parameter A. Moreover, this CFT can be put 
into a form of the harmonic fluid description for Luttinger liq- 
uids, with the Haldane controlling parameter identified with 
the statistics parameter (of quasi-particle excitations). Thus 
we propose to use the latter to characterize the fixed points 
of I-d Luttinger liquids. Such a characterization is further 
shown to be valid for generalized ideal gas of particles with 
mutual statistics in momentum space and for non-ideal gas 
with Luttinger-type interactions: In either case, the low tem- 
perature behavior is controlled by an effective statistics vary- 
ing in a fixed-point line. 



I. INTRODUCTION 

It is well-known that the Landau theory of Fermi liq- 
uids fails to describe most of one-dimensional (1-d) in- 
teracting many-body systems. To provide a substitute, 
Haldane proposed, years ago, the concept of the Lut- 
tinger liquid , defined by a set of low- lying excitations 
and critical exponents of the asymptotic correlation func- 
tions. Like Fermi liquids, there is a (pseudo-)Fcrmi sur- 
face for the quasiparticle-like excitations in Luttinger liq- 
uids, so that the classification of low-lying excitations is 
similar to that in Fermi liquids. However, the exponents 
of the asymptotic correlation functions (at low tempera- 
ture) are distinct from those for Fermi liquid theory. For 
one-component systems (without internal degrees of free- 
dom), the low-energy or low-temperature behavior of a 
Luttinger liquid is controlled by a single parameter, the 
Haldane controlling parameter. It controls not only all 
exponents, but also the velocity ratios between different 
types of elementary excitations. The Fermi liquid theory 
is a special case of the Luttinger liquids with the Haldane 
parameter A = 1. 

In recent years, the failure of Landau's theory of Fermi 
liquids to describe several newly discovered strongly cor- 
related electron systems have revived the interests in the 
theory of Luttinger liquids. Among other questions, com- 
pared with Fermi liquids, one would like very much to 



know the answer to the following questions: 

• What is the physical meaning of Haldane's con- 
trolling parameter? Or more precisely, how to use 
physical properties of low-lying excitations to char- 
acterize the concept of Luttinger liquids? 

• Does Haldane's theory of Luttinger liquids possess 
universality in one dimension, just like Landau's 
theory of Fermi liquids in three dimensions? Or 
equivalently, in terms of modern language of renor- 
malization group, does the Luttinger liquids de- 
scribe the infrared (or low-energy) fixed points in 
1-d systems? 

• In what directions could one expect to go for gen- 
eralizing the concept of Luttinger liquids to higher 
than one dimensions? 

In short, a characterization of 1-d Luttinger liquids, other 
than using a bunch of excitations and exponents, is in 
demand for gaining more insights and looking for possible 
generalization. 

To achieve this, let us recall what motivated Landau's 
concept of Fermi liquids, which is known to describe an 
infrared fixed point (or a universality class) of interacting 
electron systems. The basic idea behind it is based on 
the following organizing principle for interacting many- 
body systems: At low temperature, the low-lying excited 
states of an interacting many-body system above a stable 
ground state can be viewed as consisting of weakly cou- 
pled elementary excitations. Here "weakly coupled" only 
means that the total energy can be written as a sum of 
single-particle (dressed) energies, while the dispersion of 
the dressed energy may well depend on the total parti- 
cle number, a signal of remnant interactions between the 
quasiparticles. According to Landau, the ground state 
and the low-lying excited states of a Fermi liquid are ap- 
proximately, to a good accuracy at sufficiently low tem- 
perature, described by those of an ideal Fermi gas with 
dressed energy for the quasiparticles. 

We note the significant role played by the ideal Fermi 
gas distribution (with dressed energy) in this descrip- 
tion of Fermi liquids. Actually it is the ideal Fermi gas 
that gives a characterization to the Fermi liquid fixed 
point, and a meaning to the universality of the concept 
of Fermi liquids. This inspires us to try to give a char- 
acterization of the 1-d Luttinger liquids along a similar 
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line of thoughts, namely using a properly generalized con- 
cept of exclusion statistics, of which a special case is the 
usual Fermi statistics. Because the concept of quantum 
statistics in statistical mechanics is independent of the 
dimensionality of a system, a characterization of 1-d in- 
frared fixed points using statistics, if successful, would 
shed light on how to generalize to higher dimensional 
non- Fermi liquids. 

Fortunately, a generalization as such has been available 
recently, under the name of fractional exclusion statis- 
tics (FES). It is based on a new combinatoric rule for 
the many-body state counting ]^||, which is essentially 
an abstraction and generalization of Yang- Yang's state 
counting in 1-d soluble many-body models. FES 
has been shown to be applicable to elementary excita- 
tions in a number of exactly solvable models for strongly 
correlated systems [||,||J|-j^] , anyons in the lowest Landau 
level [^0|,p[ , and quasiparticle excitations in the fractional 
quantum Hall effect [|||J^,|ll[. The thermodynamics of 
the so-called generalized ideal gas (GIG) associated with 
FES have been studied || in a general framework. 

Inspired by these results, the thoughts along the lines 
indicated in the above paragraphs have led two of present 
authors jl2| to propose that at least for some strongly 
correlated systems or non-Fermi liquids, their low- energy 
or low-temperature fixed point may be described by a GIG 
associated with FES, similar to the way that of the Fermi 
liquid fixed point by the ideal Fermi gas (jl3|. As a tes- 
timony to this proposition, a sketchy proof was given in 
that short letter jl2| that the low-T critical properties 
of the 1-d Luttinger liquids are exactly reproduced by 
those of ID ideal excluson gas(IEG), if one identifies the 
Haldane parameter of the former with the statistics pa- 
rameter A of the latter. (We call the particles obeying 
the FES without mutual statistics exclusons). Threfore 
IEG can be used to describe the fixed points of the Lut- 
tinger liquids. In this paper, we will present our results 
obtained in |l^] in details, much of which was not pub- 
lished before. 

A main tool we use in this study is bosonization of the 
1-d excluson systems at low T, a la Tomonaga [|l4| and 
Mattis and Lieb jl5| . To bosonize an IEG system is a lit- 
tle bit tricky, because at low temperature the linearized 
dispersion of dressed energy versus pseudo-momentum 
has different slope outside and inside the pseudo-Fermi 
sea: There is 'refraction' at both pseudo-Fermi points. 
In spite of this, we still manage to construct well-defined 
density fluctuation operators that obey the U{1) current 
algebras and physically describe free phonons. Then, the 
Tomonaga-Mattis-Lieb bosonization applies, resulting a 
bosonizcd effective field theory, in agreement with Hal- 
dane's harmonic fluid description of the Luttinger liquid 
JT(| . Then the asymptotic correlation functions and their 
exponents can be systematically calculated. In this way, 
the critical properties of IEG reproduce those of the Lut- 
tinger liquids. 



An important consequence of our bosonization is that 
the low energy behavior of IEG is controlled by an orb- 
ifold conformal field theory (CFT) with central charge 
c = 1 and compactified radius |l9| R = y/T/X. This vari- 
ant of c = 1 CFT is not the ordinary c = 1 CFT compact- 
ified on a circle S 1 , rather it is compactified on an orbifold 
S 1 jZi, which is topologically an interval |l9| ]. The dif- 
ferences arise due to different selection rules for vertex 
operators, that constrain quantum numbers of possible 
quasiparticle excitations in the system. In the usual lit- 
erature this difference quite often is overlooked. Only 
within the orbifold CFT the IEG with statistics param- 
eter A = 1 recovers ideal Fermi g it should be. Also 
the two classes of c = 1 CFT's have different duality 
relation; only the one in orbifold CFT reproduces the 
known particle-hole duality in IEG, A <-> 1/A, as given 
in (For the details and more elaboration, see be- 

low. We note || that a similar situation happens for 
the Calogero-Sutherland (C-S) model The low- 

energy effective field theory for the bosonic and fermionic 
C-S models belongs to, respectively, the above-mentioned 
two classes of c = 1 CFT.) 

The fact that the low-T behavior of IEG is controlled 
by a conformally invariant theory is significant, imply- 
ing that indeed IEG provides a characterization of in- 
frared fixed points, having the conformal invariance as 
required by renormalization group. We have also studied 
the effects of mutual statistics between different pseu- 
domomenta and of the Luttinger-type (density-density) 
interactions among exclusons. In either case, the low-T 
behavior is controlled by an effective statistics X e ff for 
excitations near the Fermi points, the same way as A in 
the case of IEG. In one dimension both the momentum- 
independent part of interactions and change in chemical 
potential pL are relevant perturbations |l3|j25[| , leading to 
a continuous shift in the fixed-point line parameterized 
by A. All these will be explained in details in the present 
paper. 

To make this paper self-contained, we devote the next 
two sections, Sec. II and Sec. Ill, to reviewing the Lut- 
tinger liquid theory and the GIG associated with FES, 
respectively. In Sec. IV, we discuss the low-energy be- 
havior of the IEG system and achieve its bosonization. 
In Sec. V, the generalization to the GIG with mutual 
statistics as well as the non-ideal gas with FES are pro- 
vided. The last section is dedicated to conclusions and 
discussions. 

II. LUTTINGER LIQUID 

The Luttinger liquid, which describes a very large class 
of one-dimensional interacting many-body systems, is in- 
troduced because of the infrared divergence of certain 
vertices in the Fermi liquid description of the 1-d systems. 
Some pioneering works have been done in the Luttinger 
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model before the Luttinger liquid concept pq , [l5| . The 
model has been exactly solved by using the bosonization 
technique [Gl| . Haldane Q re-solved the model with the 
following important observations: 

(i) Besides a linearized spectrum of non-zero mode exci- 
tations, i.e., the density fluctuations (sound waves), there 
are two kinds of zero mode excitations, single-particle 
excitations by adding extra particles to the system and 
persistent currents by making Galileo boosts. 

(ii) There is a fundamental relation among the velocities 
of these three types of excitations 



(2.1) 



V s = y'VN Vj, 



where v s is the sound velocity, vj the current velocity and 
vn a velocity related to the change in particle number. 
The velocity ratios define a controlling parameter, e~ 2ip , 
by 



VN 



-2 V 



Vj 



v s e 



2<p 



(2.2) 



(iii) The above defined controlling parameter mea- 
sures the essential renormalized coupling constant, and 
is the unique parameter that determines the exponents 
of power-law decay in the zero-temperature correlation 
functions. 

Based on these observations, Haldane defined the Lut- 
tinger liquids as 1-d systems that have similar behav- 
ior (i)-(iii) at low temperature just like the Luttinger 
model. In this way the Luttinger liquids are character- 
ized through their excitations and the exponents of the 
asymptotic correlation functions. 

To be more precise, recall that the Luttinger model 
describes a one-dimensional interacting fermion system 
with the Hamiltonian 



H = J dx\\7i}\ 2 + i J J dxdyV(. 



x-y)p{x)p{y). (2.3) 



In the low energy limit, the Hamiltonian (2J5) can be 
bosonized as 

H = v s Y J \d\b\b q + \{-k/L){v n M 2 + vj J 2 ), (2.4) 



where b q are the standard boson annihilation operators, 
and M and J the operators corresponding to adding 
extra particle and boosting persistent currents, whose 
eigenvalues obey the following selection rule, 



(-1) J = (-1) 



M 



(2.5) 



The total momentum of the system also has a bosonized 
form 



P= [k F + (7r/)M]J + J2qblb q , 



(2.6) 



Eqs.(pT 2.2 J2T-2.6) turned to be universally valid for 
the description of the low-energy properties of gapless in- 
teracting one-dimensional spinless fermion systems even 
for those not exactly soluble with a conserved current 
J. This universality class is named as the Luttinger liq- 
uid by Haldane The Luttinger liquid has a model- 
independent representation, namely the harmonic fluid 
description [ fl6| |, which is convenient for calculating the 
correlation functions. The results of the harmonic fluid 
representation are listed in the Appendix, for later use to 
be compared with our bosonization theory of the IEG. 

The Haldane theory of Luttinger liquids is based on 
the significant observation that the low-T behavior of the 
Luttinger model is universal. Naturally arises the ques- 
tion: Why is it so? In this paper we intend to answer 
this question by pointing out a profound coincidence of 
the low-T behavior of the Luttinger model and that of 
ideal excluson gas (IEG), i.e., ideal gas of particles obey- 
ing fractional exclusion statistics: The universality of the 
former is due to that of the latter. 



III. GENERALIZED IDEAL GAS 

In quantum mechanics, there are two ways to define 
the statistics of particles. One is in terms of the sym- 
metry of the many-body wave function under particle 
exchange. The other is based on the state counting. 
Here we are interested in the latter definition. As is 
well-known, bosons and fermions have different count- 
ings for many-body states, or different statistical weights 
W: The number of quantum states of N particles oc- 
cupying a group of G states is, for bosons and fermions 
respectively, given by 

(G + N-lV G 1 

w »= nhg-i)i ' or W '=ww^m- (3 - 1} 

A simple interpolation between bosons and fermions is 
given by 



W = 



[G+(N-1)(1-X)}\ 
Nl [G - XN - (1 - A)]! 



(3.2) 



with A = corresponding to bosons and A = 1 to 
fermions. The physical meaning of this equation is the 
following: By assumption, the statistical weight remains 
to be a single combinatoric number, so one can count 
the states by thinking of the particles effectively either 
as bosons or as fermions, with the effective number of 
available single-particle states being linearly dependent 
on the particle number. 



G% = G-\(N-1), or 



G$ = G-(1-A)(JV-1). 

(3-3) 



with kp being the Fermi momentum. 
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Obviously, for genuine bosons (or fermions), G^ff (or 

) is independent of the particle number. In all other 
cases, either of the two G e ff is linearly dependent on the 
particle number. This is the defining feature of the FES. 
The statistics parameter A tells us, on the average, how 
many single-particle states that a particle can exclude 
others to occupy. A proper understanding of this has 
been discussed in J^?]]. Thus, the expression (^2) for 
the statistical weight, W, formulates a generalized Pauli 
exclusion principle, as first recognized by Haldane 

It is easy to generalize this state counting to more than 
one species, labeled by the index i: 



W = 



n 



i (Ni)\ [Gi-l-^XiiiNj-Sij)]] 



(3.4) 



Here Gi is the number of states when the system consists 
of only a single particle of species i. By definition, the 
diagonal Xu is the "self-exclusion" statistics of species 
i, while the non-diagonal Ay (for i =/= j) is the mutual- 
exclusion statistics. Note that Ay, which Haldane ||] 
called statistical interactions, may be asymmetric in i and 
j. The interpretation is similar to that of the one-species 
case: The number of available single-particle states for 
species i, in the presence of other particles, is again lin- 
early dependent on particle numbers of all species: 

Geff.i = - E - Sji), 



or 

G^ = G^ l + Nl -l. 



(3.5) 



The definition (3.2) or (pT4|) starts with a postulated 



form for the statistical weight, and thus is more direct 
and convenient for the purpose of formulating quantum 
statistical mechanics. One of us |3| has first formulated 
the quantum statistical mechanics by proposing the no- 
tion of generalized ideal gas(GIG): A GIG satisfies the 
following two conditions: (i) The total energy (eigen- 
value) is always of the form of a simple sum, in which 
the i-th term is linear in the particle number N^. 



E 



(3.6) 



Z=J2 Wm}) exp{]T NiQM - e°)/T} . (3.7) 
{Ni} i 

As usual, we expect that for very large Ni, the summation 
has a very sharp peak around the set of most-probable (or 
mean) particle numbers {Ni}. Using the Stirling formula, 
introducing the average "occupation number per state" 
defined by ni = Ni/Gi, and maximizing 

^- [lnW + 2 GiniOxi-e?)/T] =0, (3.8) 

1 i 

one obtains the equations that determine the most- 
probable distribution of ni 



^{SijWj +gij)nj = 1 , 



(3.9) 



with gij = XijGj/Gi, and Wi being determined by the 
functional equations 



= e (^-«)/T 



(3.10) 



The thermodynamic potential Q = —TXvlZ and the 
entropy S are then given by 



n = -PV = -Tj2Gi\og 1 ' 

i 

= -tY j g 1 H\ + w~ 1 



j 9ij n j 



(3.11) 



S : 



- In 



T 



} 



1 - 

E Gi{nfi^- + ln(l + w-')}. (3.12) 



Other thermodynamic functions follow straightforwardly. 
As usual, one can easily verify that the fluctuations, 
(Ni 2 — Ni )/Ni , of the occupation numbers are negli- 
gible, which justifies the validity of the above approach. 



with e( identified as the energy of a particle of species 



i; (ii) The state-counting ( |3.4D for statistical weight W 
is applicable. When there are no statistical interactions 
(i.e., Xij — for i ^ j), we have the usual ideal gas, 
which we call as IEG. 

With the assumptions ( |3.6|) and ( ft.4| ), the thermody- 
namics of a GIG can be worked out by the usual tech- 
niques in statistical mechanics. Consider a grand canon- 
ical ensemble at temperature T and with chemical po- 
tential [ii for species i, whose partition function is given 
by 



IV. BOSONIZATION OF 1-D IDEAL EXCLUSON 

GAS 

Let us first consider the simplest case, the 1-d IEG 
without internal degrees of freedom. We expect to ob- 
tain a continuous interpolation between the usual ideal 
Bose and ideal Fermi gas. Moreover, we want to show 
that the low-energy behavior of the IEG reproduces that 
of the Luttinger liquid and, therefore, provides a better 
characterization of the infrared fixed points associated 
with the Luttinger liquid. 
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A. Ideal Excluson Gas 

Consider a GIG of N particles on a ring with size L. 
Single-particle states are labeled by pseudo- momenta ki. 
The total energy and momentum are given by 



(4.1) 



According to ( |3.5| ), in the thermodynamic limit the hole 
density p a (k,T), (or the density of available single- 
particle states) is linearly dependent on the particle den- 
sity, p(k,T). By definition, the statistics interaction ma- 
trix is given by 



\{ki , kj ) 



Ap a (/Ci) 



(4.2) 



Or in the thermodynamic limit, one has 

X(k,k') = -Sp a (k)/Sp(k'). (4.3) 

The system is called an IEG of statistics A (with no mu- 
tual statistics between different momenta) , if 



X(ki, kj) = XS(ki - kj), 



(4.4) 



or ( |3.5|) reads 

p(kj) = J- + 1(1 - A) £ 5{k 3 ~ h)p(ki)Ak, (4.5) 



which, in the thermodynamic limit, can be simply written 
as 1 



p a (k,T) + Xp(k,T) = p (k,T), 



(4.6) 



where po{k) = l/2n is the bare density of single-particle 
states. Thus, A = 1 corresponds to fermions, and A = 
to bosons. The thermodynamic potential, now reads, in 
terms of (3.11) 



rp />00 

fl = -— dklnil + wik^)- 1 ), 

27T 7-oo 



(4.7) 



with the function w(k,T) = p a (k)/p(k) satisfying an al- 
gebraic equation, 

w(k,T) x [l + w(k,T)] 1 - x = e ( fc2 -A0/T_ ( 48 ) 

Firstly, we consider the ground state, in which the 
particles are distributed in a finite and origin-symmetric 
interval in the pseudo-momentum space. The (pseudo- 
) Fermi momentum is defined by 



kp — /i 



(4.9) 



and its value is fixed by the average particle density do 
N /L in the ground state, 



k F 



dkp(k) = c?o- 



(4.10) 



Because holes are absent in the ground state, the particle 
density in the ground state is easily obtained from (O) , 



P(k) = { 2^A' f ° r |fc| < kF; 

{ 0, for \k\ > k F . 

Hence, one has 

kp = nXdo, p — (7rAdo) 2 - 



(4.11) 



(4.12) 



Then the ground state energy and momentum are given 
by 



Eo = 
L 

Pn = 



dk P {k)k 2 



i 



7T 2 A 2 fS, 



dkp(k)k = 0. 



(4.13) 



k f 



Now let us examine possible excitations in an IEG. 
First there are density fluctuations due to particle-hole 
excitations, i.e., sound waves with velocity (see the next 
subsection) 



v s = vp — 2kp. 



(4.14) 



Besides, by adding extra M particles to the ground state, 
one can create particle excitations, and by Galileo boosts 
a persistent current. We observe that the velocities of 
these three classes of elementary excitations in IEG also 



satisfy the fundamental relation (2.1). Indeed, shifting 
A^o to N = No + M, the change in the ground state 
energy is 

S M E = \ 2 X 2 {N/Lf - ^ 2 X 2 {N /Lf 

= tt 2 (P M + n{Xk F )M 2 + 0(M 3 /L 3 ), (4.15) 

while a persistent current, created by the boost of the 
Fermi sea k — > k + irJ/L, leads to the energy shift 

/kp-\-7rJ/ L rkF 
dkp(k)k 2 - / dkp(k)k 2 
-kp+TrJ/L J —kp 



ir(k F /X)J 2 



(4.16) 



Therefore the total change in energy, due to charge and 
current excitations, is 



5£ - fj,M = —v F {XM 2 + A" 1 J 2 ). 
21j 



(4.17) 



The total momentum change due to the current excita- 
tions is 



SP = 22-r = n ( d o + -r)J- 



(4.18) 
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If we denote the variation in free energy as SFq = 8Eq — 
fxM, and identify A as the controlling parameter e~ 2ip in 
the Luttinger liquid theory, ( 4.17[) j ust recuperat es the 
zero-mode contributions in (|2.4|). Comparing (4.17) 
with (E^J) we identify the velocities vn and vj to be 



VN — W_fA, Vj = Vp/X, 



(4.19) 



Then we see the velocity relation (p.l|), i.e., v s = y/vjpvj, 
that Haldane used to characterize the Luttinger liquids, 



is satisfied in IEG. The selection rule {2.1) also holds for 
the IEG, since the system should correspond to the ideal 
Fermi gas if A = 1. 

Encouraged by this relationship between the IEG and 
Luttinger liquids, we want to calculate the critical expo- 
nents of IEG to see whether they reproduce those of the 
Luttinger liquids. This motivates to develop a bosoniza- 
tion for the density fluctuations in IEG. 



B. Low Energy Limit and Bosonization 



Following Yang and Yang jl 18 1, we introduce the 
dressed energy e(k,T) by writing 



,{k,T) 



_ e(k,T)/T 



(4.20) 



The point is that the grand partition function Zc, cor- 



responding to the thermodynamic potential (4.7), is of 
the form of that for an ideal system of fermions with a 
complicated, T-dependent energy dispersion given by the 
dressed energy: 



z G =na+e- e(w ). 



(4.21) 



However, this fermion representation is not very useful, 
because of the implicit T-dependence of the dressed en- 
ergy. To simplify, we consider the low-T limit. By using 
the dressed energy, (4.8) reads 



e(k, T) = k 2 - fx - T(l - A) ln(l + e - e(k ' T)/T ). (4.22) 

Because there is no singularity in e(k,T) at T — 0, the 
zero temperature dressed energy is given by 



<k) = 



{k 2 -k 2 F )/X, \k\ <k F , 
k 2 - kp, \k\ > k F . 



Denote 

e(k,T) = e(k)+e(k,T), 

where 

e(fc,0) = 0. 
In the low-T limit, one has 



(4.23) 



(4.24) 



(4.25) 



e(k T) = J - ( A_1 - i) 7 ^ 1 + e- |£(fc)l/T ), \k\ < k F , 

[ ' ' \ (fc 2 - M )-(l-A)7Tn(l + e-l £ ( fe )l/ T ), |fc| > k F , 

(4.26) 



Hence, 



m = 



(1 - A- 1 )Tln(l + e-l e «l/ T ), |fc| < k F , 
(A-l)Tln(l + e-l £ ( fe )l/ T ), |fc| > k F . 

(4.27) 



For low energies, one can consider only the excitations 
around the Fermi surface, 



Q(T) 



/ dklnll 



-e(fc)/T 



1 

2^ 
1 

2^ 
T 



dk- 



e(k,T) 



— OO 

k f 



— kp 



dke(k) 



oe(k)/T 

T 



7rA 



dfcln(l + e 



k F +S 



dk ln(l 



kp — 8 
(*)\/T\ 



-\<m/T\ 



(4.28) 



where the first term on the right hand side of the last 
equality is recognized as ft(0)/L. The cut-off S is of order 
0(T/v s ) (actually, a few times of T/v s ). Mathematically, 
we take the limit of T — > followed by S — > 0. Using the 
integral formula 



dxln(l + e~ x ) 



12' 



we have the low-T thermodynamic potential 
Q(T) n(0) irT 2 



L 



L 



6v, 



(4.29) 



which implies that the theory is cut-off independent at 
low temperature. Notice that F = SI — /iiV '. Because 
we only consider the particle-hole excitations near the 
Fermi surface contribute to thermal excitations, N(T) — 
N(0) = 0, which can be checked by an explicit calculation 
in terms of the definition of p(k, T). Thus, we have 



F(T) F(0) Q(T) S7(0) 



L 



L 



L 



irT 2 
6v„ 



(4.30) 



This means that the low energy behavior of the IEG is 
controlled by a c = 1 CFT. This result can be verified by 
a finite-size scaling in the spatial direction, 



F L (0) F{0) 
L 



L 



(4.31) 



where Fl(0) is the zero temperature free energy for a 
system with size L. (For details, see Q.) 

The above relation agrees with the finite-size scaling of 
a conformally invariant system with central charge c = 



6 



1. So we want to see whether the low-energy effective 
theory of the IEG is really a CFT. Let us start with the 
grand partition function ( 1.21] ). At low temperature, the 
solution J4.26D leads to e(k, T) = 0{Te~ ^ /T ), so one can 
simply replace e(fc,T) with e(fc) in the grand partition 
function: 



)• 



(4.32) 



Note that the dressed energy with k outside the Fermi 
points ±kp has a slope different from that with k inside 
±k F . The former is ^7 and the latter It is nec- 

essary to keep this in mind for writing down the correct 
ground state wave functions of the excluson system. Now 
that the dispersion e(k) is T- independent, the grand par- 
tition function in the low-T limit can be expressed in a 
fcrmionic representation as 



Tre 



where the effective Hamiltonian is given by 
Has = e(k) c[ck, 

k 



(4.33) 



(4.34) 



where ct are fermionic creation operators. We also see 
that e(kp) = 0, which can be used to define the Fermi 
momentum. 

Physically, it is the phonon excitations that dominate 
the low-energy behavior of the system. In the low-T 
limit, it is enough to consider the density fluctuations 
only near the Fermi points, k ~ ±kp, where the left- and 
right-moving sectors are separable and decoupled: 



(4.35) 



Besides this, another important simplification for excita- 
tions near Fermi points in the low-T limit is that their 
energy, H±, has a linearized dispersion: 



e ± (k) 



±vp(k^fkp), \k\>kp, 
±vp(k =p kp)/ X, \k\ < k F . 



(4.36) 



We note the 'refractions' at k = ±kp, which implies to 
create a particle with pseudo-momentum k and to create 
a hole with kl cost different energies, even if \k — kp\ — 
I As' — kp\. The reason for this is that k is not the actual 
momentum carried by ct, as we will see soon. 

The key thing for bosonization is to construct a density 
fluctuation operator. Taking into account the different 
slopes for dressed energy inside and outside the Fermi 
points, the density fluctuation operator at k ~ kp is 
constructed as follows: 



P q +) = 12 : c l+g Ck : + 12 : c l+\g c k 

k>kp- k<kp- — Xq 



12 



kF — Xq<k<kF 



c \~k F Ck ■ 



(4.37) 



for q > 0. A similar density operator p q can also be 
defined at k ~ — fcp, 



P [ q- ] = 



12 )C k-q C k 

k<-k F 

12 < 

kF-\-Xq>k> — kp- 



■ 12 A,,.'/. 

k> — kp+Xq 



k+k F C k ■ 



(4.38) 



To define the normal ordering we write, e.g. 

_ J c k , k > kp, 
I d\, k < k F , 



(4.39) 



where dl is understood as a creation operator of a hole. 
Then normal ordering is done as usual: putting the anni- 
hilation operators to the right of the creation ones. Hence 
we have, e.g., 



Hq 



12 :C Uq C k 
k>k F 



12 



-Xq<k<kp- 



- k-k F . ,^k 
— ^ +q 



12 '■ d k+\q d { '■ 
k<kf —Xq 



(4.40) 



Within the Tomonaga approximation p9[ , in which 
commutators are taken to be their ground-state expecta- 
tion value, we obtain 



J2 (0\c k +x q ct +Xql \0) 

k p — Xq<k<k f 

= S q,q' 12 1 = 2^ qS 1'4' 

kp —Xq<k<kp 



(4.41) 



Also, the commutators between H e g and p q are 

[H ± ,pW] « (0\[H ± ,pW}\0) = ±v F qp[ ± \ (4.42) 



( 1.41] ) and ( 4.42Q describe 1-d free phonons with the 
sound velocity v s — vp (so we have proved ( 4.14 )). In- 
troducing normalized boson annihilation operators 



b q = yj2n/qLp q - 



(+) 



I qLp^ (4.43) 



and adding back the zero mode contributions, the 



bosonized Hamiltonian satisfying (4.41) is given by 

1 
A 



H B = v s {J2 l(b\b q + b\b q ) + l^-[\M 2 + I J 2 ]}, 

g>0 



2L' 



(4.44) 



which agrees with the bosonized Hamiltonian (2.4) in the 
Luttinger liquid theory. 

In passing, we make a comment on linearization of the 
dressed energy dispersion. When we did this, we changed 
the ground state energy, because we assumed that for 
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all k the spectrum is linear in k. However, we changed 
neither the ground state wave function, nor the low-T 
physics. On the other hand, the linearized spectrum was 
valid only for phonon excitations, it has nothing to do 
with the zero-mode excitations. So, after the linearized 
phonon part of the Hamiltonian is bosonized, we had to 
add back the zero- mode excitations. 

The construction of the bosonized momentum opera- 
tor is a bit more tricky, because c k does not carry a mo- 



mentum k. Each term in (4.37) should carry the same 
momentum q, therefore the fermion created by ct carries 
a dressed momentum p, which is related to k by 



p(k) 




k > kp, 
\k\ < k F , 
k < —kp. 



(4.45) 



In terms of this variable, the linearized dressed energy 
e(p) is of a simple form: e±(p) — ±v s (p =F Pf), with 
Pf = kp/X. The bosonized total momentum operator, 
corresponding to the fermionized P = ^2 k p(k) c\,Ck, is 

^ = E - ~ b \ l i) + ^ + M / L ) J - ( 4 - 46 ) 

q>0 

We see that the fundamental velocity relation, the 
bosonized Hamiltonian and momentum, and the selec- 
tion rule of the quantum numbers in the Luttinger liquid 
theory can all be reproduced in IEG if we identify 



A 



-2tp 



(4.47) 



To say that IEG can be used to characterize the renor- 
malization group fixed points of Luttinger liquids, we still 
need to check the conformal invariance of the bosonized 
theory of IEG, and to verify the critical properties of IEG 
reproduce those of the Luttinger liquids. 



C. Effective Field Theory and Conformal Invariance 

To check conformal invariance, we need to rewrite 



the above bosonized effective Halmitonian (4.44) into a 
form of field theory in coordinate space. Employing the 
Fourier transformation, the density operator can be writ- 
ten as 

p{x) = p R (x) +p L (x), 



<J>0 



where Mr,l are given by M = Mr + Ml and b q = b- q 
for q > 0. 



The boson field </>(x), which is conjugated to p{x) and 
satisfies 



[<j)(x),p(x')} = iS(x-x'), 

is given by 

<t>(x) = (f> R (x) + (f>L{x), 

9r{x) = ~z~ 



(4.49) 



g>0 



4>a . t^Jrx 



<p L (x) = y + ^jr +*E y ^ iqx ~ b « - e<9 ^)> ( 4 - 5 °) 

with J = Jr + Jl- We have to assign the quantum 
numbers such that there are only two independent each 
other in Mr^l and Jr,l- A consistent choice is 



Then, 



Mr = J R , M l = -J L . 



J = Jr + Jl, M = Jr-J l . 



(4.51) 



(4.52) 



Here (j)Q is an angular vari able conjugated to M: 
[<f>o,M] = i. The Hamiltonian (|]J) becomes 

H = \j Q dx[nv N p(x) 2 + vj/n (d x cj)(x))% (4.53) 
or by a field rescaling, 



H=£-\ dx [n(xf + (d x X(x))% (4.54) 
^ Jo 



where 



IL(x) = 7rX^ 2 p(x), X(x) = \- x / 2 (j){x). (4.55) 

With X(x,t) = e tHt X (x)e~' iHt , the Lagrangian den- 
sity reads 



C= - L d a X{x,t)d a X(x,t). 

Z7T 



(4.56) 



This is the Lagrangian density of a free scalar field theory 
in 1 + 1-dimensions. Writing the corresponding operators 
as the functionals of the scalar field X(x,t), all correla- 
tion functions can be obtained by using the propagators 
of XR(x,t) and XL{x,t), 

(X R (x, t)X R {0, 0)) = ~\n.(x - v s t), 

(X L [x,t)X L (0,Q)) = -~la(x + v a t). (4.57) 

The statistics of an operator in the theory can also be 
inferred by the commutators of the scalar fields, 



[X RtL (x),X RtL (x')] = ±—9(x - x'). (4.58) 
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We recognize that C (4.56) is the Lagrangian of a c = 1 
CFT consistent with the finite-size scaling (4.30). 



Alternatively, it is easy to check that the theory is in- 
variant under the conformal transformations generated 
by a set of the Virasoro generators 



- y 

2 ^ 



x-n — ni lJL n ; 



Lm, — 



1 oc 



"■n—m^n i 



(4.59) 



where the oscillators a r , 



m}/ 2 b q and cL 



= to 1 / 2 ^ for 

to = qL/2n > being integers. a = {tt/2L) 1 ^[J\- 1 / 2 - 
MA 1 / 2 ] and a = (tt/2L) 1 / 2 [J\~ 1 / 2 + MA 1 / 2 ]. The gen- 
erators obey the Virasoro algebra with the central charge 
c=l, 



(4.60) 



with Iff =L m +L m . 

Since 0o is an angular variable, there is a hidden in- 
variance in the theory under <j) — > </> + 27r. The field X is 
thus said to be "compactified" on a circle, with a radius 
that is determined by the exclusion statistics |],^l| : 



X~X + 2irR, R 2 = l/\. 



(4.61) 



Noting the selection rule (2.5), the Hamiltonian has a 
duality 



A^l/A, M<->J, 



(4.62) 



which has referred to the particle- hole duality [@J| . Us- 
ing the CFT terminology, this duality is represented as 
the duality of the compactified radii, 



R <-> 1/R. 



(4.63) 



We note that this is different from the duality relation 
R <-> 2 /R, in the usual c = 1 CFT [|l9| compactifed on a 
circle. Actually, according to the standard terminology 
in C FT JlSj ] , our selection rule ( |2.5| ) and duality rela- 
tion (4.63) make what we obtained above a c = 1 CFT 
compactified on an orbifold S 1 /Z2, i.e., a circle folded 
by a reflection about a diameter, which topologically is 
a semi-circle or an interval. This difference can also be 
seen from the grand partition function: Using the iden- 
tification between H e g and Tp ot , i.e., H c s = v s LQ 0t , the 
grand partition function of IEG (in the low-T limit) can 
be rewritten as 



Tr H [q Lo q Lo 



(4.64) 



where 



e iv ' T with t = i/3 = i/T. Thus, the selec- 
tion rule (2.5) severely constrain the allowed values for 
the eigenvalues of Lq and Lq. It makes the CFT we ob- 
tained have an unusual spectrum and duality relation, 



corresponding to the c = 1 orbifold CFT |19| . In next 
subsection we will see that because of the difference in 
the selection rules, the statistics of the allowed charge-1 
operators in the two classes of CFT's are not the same. 

We note that a similar situation happens for the 
CFT that describes the low-T behavior of the Calogero- 
Sutherland (C-S) model 0Q. This model has two 
different versions, with the long-range interactions be- 
ing among bosons or among fermions, respectively. At 
low temperature, the two versions have different selection 
rules for the zero-mode quantum numbers, thus leading 
to different CFT's: The low-T CFT for the bosonic C-S 
model is the usual c = 1 CFT compactified on a cir- 
cle, which has been studied extensively in the literatures 
foj-flJ); while for the fermionic C-S model the low-T 
limit gives rise to the c = 1 orbifold CFT. This is be- 
cause the selection rule for zero modes severely constrains 
the spectrum of the system, i.e., possible quantum num- 
bers of the allowed excitations. (For details, see ref. 
Thus, only the fermionic (not the bosonic) C-S model 
respects a duality relation A <-> 1/A that coincides with 
the particle-hole duality in IEG 



D. Correlation Functions 

The CFT description of the IEG offers a better under- 
standing for the space of quantum states in the theory. 
States V[X] |0) or operators V[X] are allowed only if they 
respect the invariance (4.61), 



V[X + 2irR] = V[X], 



(4.65) 



with a given boundary condition restriction. Here, a 
Fermi or a Bose operator obeys the periodic boundary 
condition (PBC). So quantum numbers of quasiparticles 
are strongly constrained, in particular by the selection 
rule for zero-mode quantum numbers. For example, the 
primary fields obeying the PBC in the CFT are given by 

<Pau(x) ~ /( J, X°) : e i(M*"+J/#' a )x*(*) 

t(M\^ 2 -J/\^ 2 )X L (x) . 



x e 



f(J,X»)=e 



0) = jj(\ 1/2 -\- 1/2 )x% -iJ(X^ 2 -\-^ 2 )X° 



(4.66) 



where the prefactor f(J,X°) makes the fields satisfy the 
PBC, M and J eigenvalues of the number and current op- 
erators, and X° = irMx/L. The field carries the charge 
M and current J. The conformal dimensions of the fields 
are 



h = ^[(MA 1 / 2 + J/A 1 / 2 ) 2 + (MA 1 / 2 - J/A 1 / 2 ) 2 ] 



= M 2 A + J 2 A" 



(4.67) 



The statistics of the field can be calculated by using 
(4.58) and the statistics factors are 
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exp{» JpfA 1 / 2 + J/X 1 ' 2 ) 2 - {MX 1 ' 2 - J/X 1 / 2 ) 2 }} 

(4.68) 



-1 



,AfJ 



Consider the charge-1 primary fields, with M = 1. 
Therefore, they can only be fermions since J = odd due 
to the selection rule. The general charge-1 fermion op- 
erator is a linear combination of the charge-1 primary 
fields. A careful construction of the allowed fermion field 
with unit charge leads to 



anyon field needn't to obey the PBC. So in the theory, 
there may be allowed anyonic excitations. A charge-1 
anyonic (or exclusonic) operator is a primary field that 
does not obey the PBC, 

*{{x) =: ^ F ( x ) e ^ 1/2 -*- 1,2 )(Xn{x)-x L(x) ) . (4 ?4) 
The anyon commutation relation is easy to check: 



$>{{x)*[{x') - e inXssn(x - x 'HUx')^{(x) = 0, 



oo 
m— — oo 



a i{\ 1 ' 2 + (2m+l)/\ 1 ' 2 )X R (x_) . 



for i^i' 
(4.75) 



: e 



(\ 1/2 -(2m+l)/\ 1/2 )X L (x + ) . 



In other words, the anyon field carries a fractional cur- 
(4 69) ren t- O r by the M <-> J-duality, the anyon with integer 



where the prefactor / has been suppressed and the hcr- 
mitian, constant-valued operators O m satisfy |30f| 



J carries a fractional charge. The correlation function of 
the single-anyon reads 



[O m , O m >] = in(m - m'). 



(4.70) 



G(x,t;X) = (^{(x,t)^ x (0,0)) 



The multi-sector density operator is the linear combi- 
nation of those primary fields with M — and J =even, 



■ £ 

m— — oo 



m (m+A) 2 /A m 2 /A 
X R X L 



(4.76) 



p(x) = W F (x)^ F (x) 



This correlation function coincides with the asymptotic 
one Q in th e Cal ogero-Sutherla nd mo del. We see that 

ex U2m\X (x) - X ( )]/X 1/2 \ • (A 71^ if A = 1; consists with d^) i 

Pi L R\ ) l{ )\l i ■ ■ \ ■ there are no boson excitations (A = 0) because 

G(x,t;0) = 0; 



All the secondary fields in the CFT follow by consider- 
ing the sound wave contribution to the conformal weight 
of the fields. 

The correlation functions can easily be calculated by 
using the CFT techniques. For examples, the density- 
density and single particle correlation functions are as 
follows, 



(p(x,t)p(0,0)) 



(iii) and moreover, A > is implied since (4.77) will 
diverge at the long distance if A < 0. 

(iv) Look at m = 0. The critical exponents can be reads 
out, 



A + A -1 , 



2A. 




Thus, 



m 2 /A 



COs(27TC?o™^) 



(4.72) 



and 



G(x,t) = (tf p (x,t)9 F (0,0)) 



d 



m— — oo 



Br, 



1 



(A 1 / 2 + (2m+l)A- 1 / 2 ) 2 /4 



1 



(A 1 /2_(2m+l)A- 1 /2)2/ 4 
X L 

o i(2Tr(m+\/2)d x+fj,t) 



(4.73) 



where xr.l = x =F v s t and A m and B m regularization- 
dependent constants. 

Usually a physical quantity, e.g., a boson field, satis- 
fies the periodic boundary conditions (PBC). Hence, a 
charge-1 bosonic excitations are not allowed in the the- 
ory, because it is anti-periodic. However, as we know, an 



rjf > 7] X , if A < 1; 
T]f < 77A, if A > 1. 

(v) The multi-sector density operator for exclusons is the 
same as that of the fermion. 

The single-hole state, i.e. *t/ A |0) = * A (A -> A" 1 )^), 
with charge — 1/A alone is not allowed. The minimum 
allowed multi-hole state is given by 

^(^•••^/aWIO) 
if A = p/q is rational. One may obtain, e.g. , 

([V{ /X (x,t)r[y 1/X (0,0)] p ) ~ [G(x,t;l/X)f. 

A more interesting allowed operator is what creates q 
particle excitations accompanied by p hole excitations: 

h(x,t) = [*{(x ) t)n*\ /x (x,t)r. 

We note the similarity of this operator to Read's or- 
der parameter |31| for fractional quantum Hall fluids (in 
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bulk). Its correlation function can be calculated by using 
Wick's theorem: 

(n(x,t)n(0,0)) ~ [G(x,t; X)] q [G(x,t; 1/A)] P . (4.78) 

If the contribution from the m = sector dominates, 
then one gets 

(h{x, t)h(0, 0)) ~ (x — v s 



V. TWO EXTENSIONS 

Now we proceed to go beyond IEG. Two extensions 
will be discussed in this section: The one-component 
GIG with the mutual statistics, and the non-ideal gas 
with the Luttinger-type interactions. In either case, we 
will show that the low-temperature behavior is that of an 
IEG, controlled by a single "effective statistics" param- 
eter X e ff, whose value depends on the mutual statistics 
and the coupling constants in the interactions. 



A. Generalized ideal gas with mutual statistics 

We turn to discussing the effects of mutu al st atistics. 
Consider a GIG with the statistics matrix ( |4.3| ) in mo- 
mentum space given by 



g(k - k') = S(k - k') + $(fc - k'). 



(5.1) 



Here $(fc) = $(— k) is a smooth function. $(fc — k') 
stands for mutual statistics between particles with differ- 
ent momenta; for IEG $(fc) = (A — l)S(k). The thermo- 
dynamic properties of GIG is given by eq. (4.7 ), but now 
w(k, T) satisfies integral equation ||,pl which, in terms of 
the dressed energy (4.20), is of the form 



e(k, T) = eo(fc) +rf^§(fc- k') ln(l + e ~«W T ) , 

(5.2) 

where eo(fc) = k 2 — p. In the low-T limit, it can be proven 
by the iteration § that e(k,T) = e(k) + 0(T 2 /v s ) 1 where 
e(k) is the zero-temperature dressed energy given below. 
At T — 0, the Fermi momentum k F is determined by 



e{±k F ) = 0. 



Introduce 

{a ■ f3)[-k F , k F ] ee 

($ • a)(k; — kp, kp] = 





dk 


1 ~k F 


2^ 


rkp 


dk' 


l-k F 


2^T 



$(fc- k')a(k'). 



(5.3) 

(5.4) 
(5.5) 



Then both p(k) and e(fc) in the ground state satisfy an 
integral equations like 

a(k) = a (k) -($-a) (k; -k F ,k F ] . (5.6) 
The dressed momentum p(k) is related to p{k) by 

dp(k) = 2irp(k)dk, p(k) = -p(-k). (5.7) 

The ground state energy is given by 

E /L = (e -p)[-k F ,k F }. (5.8) 

Using the equation satisfied by p(k), it can be expressed 
by the dressed energy 



E /L = (e ■ p a )[-k Fl k F }. 



(5.9) 



The above equations are of the same form as those 
in the thermodynamic Bethe ansatz 0], hence the 
Luttinger-liquid relation [f32"[ , v s = ^/vjpuj, remains true. 
A simple proof is sketched as follows. The sound velocity 
is well-known: 



v s = de(p F )/dp F . 
The charge velocity is given by 

v N = v s z(k F y 2 , 



(5.10) 
(5.11) 



where the dressed charge z{k) p^] is given by the solution 
to the integral equation 

z(k) = l-($-z){k;-k F ,k F ]. (5.12) 

This relation can be easily derived from the definitions 



vn = Ldp/dN , 
z(k) = -6e(k)/6(i. 



(5.13) 



To create a persistent current, let us boost the Fermi sea 
by 



± k F — > ±fcf 



(5.14) 



where A = z(k F )/ Lp(k F ). Then the total energy of the 
state with the persistent current is 



E A /L = (e ■ pA)[-k F + A, k F + A] 
= (ca -po)[-k F + A,/c F + A], 



(5.15) 



where 

PA(k) = p (k) - ($ • PA )(k; -k F + A,k F + A] (5.16) 
and 

ca(*0 = eo(fc) - (* ' £A)(fc; -k F + A, k F + A]. (5.17) 

Now, using the last expression for E A and substituting 
e a hi (VA), we have 
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E A /L= ( e . po )[-k F ,k F ] + — -e'(k F ) 

{Po(&f) + (po ■ 2wF)(k F ; -k F , k F }} 

A 2 , 

-^-e'(-k F ){p (-k F ) + (pa ■ 2nF){~k F - -k F , k F }}. (5.18) 
Here F(k,k') is determined by 



P(k) 



= 1 - / dk'{$(k- k') 
Pa J-k F 

-($ • $)(fc'; -k F , k F ]($ • ($ • $)(fc'; —k F , k F ] 



(5.26) 



F(k,k>) = ±$(k,k')-±- 

Z7T Z7T 



dk"${k,k")F{k",k'). 

(5.19) 



On the other hand, we note that the equation for po(k) 
can be rewritten as 

P (k) = Po (k) - (p ■ 27rF)(fc; -k F , k F ]. (5.20) 

Thus, we have 

E A -E = LA 2 e'(k F )p(k F ) = (2tt / L)v s z{k F f . (5.21) 



This verifies vj = v s z(k F ) 2 . In view of eq. ( 4.19 ), at low 
energies, the GIG looks like an IEG with 



\ eff = z {k F y 



(5.22) 



It can be shown that it is the effective statistics ( 5.22| ) 
that controls the low-T critical properties of GIG, as A 
does for IEG. Linearization near the Fermi points and 
bosonization of the low-energy effective Hamiltonian go 
the same way as before for IEG. The only difference 
now is that the slope of the linearized dispersion for the 
dressed energy e±(k) = ±e'(k F )(k^fk F )+p = ±v s (p(k)^f 
p F ) + p, is smooth at k ~ ±k F . So bosonization is stan- 
dard and the bosonized Hamiltonian is the same as eq. 



(4.44) for IEG, only with A replaced by A e //. However, 
before going to the bosonization we need an effective 
Hamiltonian of the fermions with the dressed energy. Un- 
like the IEG, in the GIG case, e(k,T) = e{k) + 0{T 2 /v s ). 
Now, we work out the T-expansion of e(k, T) explicitly 
in the low-T limit: 



e(k,T) = e(k) + i(k,T) + 0(T 3 /v 2 s ). 
One finds that 

nT 2 



e(k,T) 



/(*). 



6e'(k F ) 

with the function / determined by 

f(k) = $(fc F - k) - ($ • /)(*; -k F , k F ] 
= $(fc F - fc) - ($ ■ $)(fc; -k F , k F ] 
+($ • ($ • <S>)(k;-k F ,k F ] + ... 



(5.23) 



(5.24) 



(5.25) 



Note that the equation that p(k) obeys can be rewritten 



Integrating (5.25) over k and comparing with (5.26), one 
has 



k F 2tt po 



and then 



dk. 



irT 2 

_ k 2tt 6e'(k F ) 



(5.27) 



(5.28) 



Substituting (5.23) into the thermodynamic potential 
(4.7), we have 



n(r) 



T 

27T 



dfcln(l + e- £ W/ T ) + 



i r°° dk 

(5.29) 



e(fc,T). 



In the low-T limit, the first term in the last equation 
gives 



n(o) 



ttT 2 



6e'(k F y 



with 



n(o) 



i 

27 



dke(k). 



and the second term is approximately given by ( 5.28 ) 
Thus, 



fi(T) n(0) ttT 2 



6v s 



(5.30) 



which proves the central charge c = 1 CFT behavior of 
the theory at the low energy. 

We may also confirm this from the finite size scaling in 
the spatial direction. To see this, we consider the discrete 
version of the equation in which the density pL{ki) obeys 



PL(h) = — -^$(fc, - kj). 



(5.31) 



Using the relation between discrete sum and integration 

N 2 j ,(JV 2 -l/2)/L 



i N * t r(N 2 -l 
L „-a,, L J JVi+l/ 



dxf(x) 



n=JVi 
-J_| 

0(1/L 3 ), 



(iVi + l/2)/i 

/'((JV 1 -l/2)/L)-/'((JV a + l/2)/L)] 



(5.32) 



one has 
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PL(k) » ($ • Ph){k; -k F ,k F ] 

Ztt 



1 



1 



24L 2 p(fc F ) 

1 1 

24L 2 p(k F ) 



d$(k - k') 
Ik* 

d<f>(k-k') 
dk' 



temperature-dependent spectrum. To derive the low- 
energy effective theory, however, one rewrites the ther- 
modynamic potential ( [5.29 ) in the low-T limit as 



J -k f 



O(T) fi(0) 



Denote 



PL = P + Pi , 



(5.33) 



(5.34) 



L 



-2Tp(k F )I(k F ,T), 



(5.39) 



where p(k) is of the order 0(1/ L°) and pi(k) the order 
0(1/L 2 ). Then, p(k) is as defined and pi{k) is deter- 
mined by 



Pi(k) 



1 



1 



24Z, 2 p(k F ) 



d$(k-k') 



dk' 



where 
I(k F ,T) 

That is, 

n(r) _ 



kF+S dk\n{l + e-^' T ) 

k F -8 

^p(p)ln(l+ e -l £ «^l/ T ). (5.40) 



Pf-6 



2tt 



feF dk CM 



24L 2 p(fc F ) 



9$(fc - fc') 



(*-Pi)(fc; (5.35) 



T 
27 



2tt 

Pf-<5 



e?pln(l + e 



-|e(fc(p))|/T 



), (5.41) 



-p F +S 



The corresponding thermodynamic potential reads 

ni(0, 4i><4» 



where p is the physical (dressed) momentum. The grand 
partition function reads 



IK* 



-/3e(fc(fc')) 



(5.42) 



dkp L (k)e (k) 



1 



-A* 

fejr 



dkp(k)eo(k) 



24L 2 p(k F ) 

dkp!(k)e (k) 



t'o( k )\-k F -e' (k)\ kf 



where k! = k for |fc| < k F — 5 and k' = p for \k\ > k F — 8. 
Now, we can have an effective Hamiltonian because e(k) 
is T- independent, 



H eS = $>(fc(fc'))4, 



Cfc'. 



(5.43) 



24L 2 ( o(fc F ) 



(5.36) 



The first term of the last equation is f2(0)/L and the rest, 
using (5.35), can be written as 



d 



24L 2 p(k F ) dk 



e (k) 



-!)($• e<>) 



+ (-l) 2 (($-$)-e )-|-...)(A : ;-A :F ,fc F 
1 9 



k—kf 



Similar to the IEG case, the low-T excitations can 
be considered by taking the linear approximation near 
the Fermi points, and after bosonization the zero- 
temperature excitations should be added back. The way 
to bosonizc the linear Hamiltonian is also similar to the 
case of IEG. Because the dressed energy is smooth at the 
Fermi points now, the bosonization is even simpler. The 
density fluctuation operators are simply given by 



2AL 2 p(k F ) dk 



e (k) + (-!)($ -e ) 



+ (-l) 2 (($.$). eo ) + ...)(fc;-A:F,^ 



(5.37) 



Pq 



k— — k.F 



k^kp 

E : Cp-qCp ■ 
k^ — kr 



Recall the equation that e(k) obeys, one has immediately, ijv^g commu t a t ors 



among p q ^ and H± are 

E 



fii(o) n(o) 



7T e'(k) kF - e'(fc)_fcj 



L 



L 



12L 2 



2np(ki. 



■TVs 

6L 2 ' 



(oIIp^pPJIo) 



(5.38) 



5g, q > 



as desired. 

Similar to the case of IEG, we also could have a fermion 
representation of the grand partition function with the 



E 



p F —q<p<PF 



PF-q<p<PF 
L 
2^ 



q5q,q> 



(5.44) 



\ c p+q c p+ q ' |0) 



(5.45) 



and 
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[H±, P ^] « (0\[H±,pW]\0) = ±«i,gpW. (5.46) 

Introducing the normalized bosonic annihilation opera- 
tors 



b q = y/2ir /qL Pq +\ b q = y/2ir/qLp( q 



(-)t 



(5.47) 



and adding back the zero-mode contributions, the 
bosonized Hamiltonian satisfying ( 5.46| ) is given by 



H B = v s {^q{b\b q + b\b q ) 

q>Q 



Z L A e ff 



(5.48) 



which agrees with the bosonized Hamiltonian (2.4) in the 



Luttinger liquid theory. We see that with A replaced by 
A e ff, the bosonized Hamiltonian for the GIG is the same 
as that for the IEG. So, all consequences we have ob- 
tained from the bosonized Hamiltonian in the IEG case 
can be applied to the GIG case. Especially, there is an 
(allowed) describing the particle excitation near 

the Fermi surface with both anyon and exclusion statis- 
tics being A e / / . In this sense, one may say that the effect 
of mutual statistics is to renormalize the statistics ma- 
trix. 

Here we remark that in IEG, $(fc, k') = (A— 1) S(k-k') 
is not smooth, so the dressed charge has a jump at fc^: 
z{k%) = 1 and z{k F ) = A" 1 for k$ = k F ± 0+ . The 
general Luttinger-liquid relation is of the form 

v N = v s [z(kj,)z(kp)}~ 1 , vj = v s z(k%)z{kp). (5.49) 



Using the Bogoliubov transformation, the Hamiltonian 
can be easily diagonalized 

H = ^io q {a\a q + a\a q ) + ^(tt/L)[v n M 2 + vjJ 2 } + S , 



q>0 



So = X> 



v s q), 



(5.52) 



q>0 



where 



a q = cosh (po b q — sinh (p^ b q , 
a q = cosh <f>o b q — sinh (p^ w q , 

and the renormalized velocities are 



(5.53) 



v s 
v N 



\(v s + U ) 2 



v N = v s e 2vo , 

J2>f>0 



vj = v s e 

with the controlling parameter tpo determined by 
vj - v N - 2V 



tanh(2(^o) 



vj + v N + 2U ' 



(5.54) 



(5.55) 



Thus, the Luttinger-liquid relation ((4. IS) survives with 
X e f f of GIG renormalized to 



A 



eff 



-2<p 



(5.56) 



Note that the new fixed point depends both on the posi- 
tion of the Fermi points and on the interaction parame- 
ters Uq and Vq, leading to "non- universal" exponents. 



B. Non-ideal Gas 



VI. DISCUSSIONS AND CONCLUSIONS 



Finally, we examine non-ideal gases, e.g., with general 
Luttinger-type density-density interactions, 

H = -ffcff + Hi , 

Hl = I ^PiiPA + PlPq) + V l(PlPl + hP\)l ( 5 ' 5 °) 
q>0 



where H c g is given by ( 5.43 ) describing a GIG, and p q 
and p q are the excluson density fluctuations near ±k F 
respectively. After bosonization, the total Hamiltonian 
remains bilinear in densities: 

H = H B +H I 

= \Y, 9[(v. + U q ){b\b q + b\b q + b q b\ + b q b\) 



q>0 



lvr, 



+V q {b\b\ + b q b q + fttftt + ~ bqbq )} + ~[v N M 2 + vjJ 2 } 
+j[Uq(M% + M 2 L ) + 2V M R M L ] - v s g (5.51) 



9>0 



In conclusion, we have shown that 1-d IEG (without 
mutual statistics) exactly reproduces the low-energy and 
low-T properties of (one-component) Luttinger liquids. 
This gives rise to the following physical picture: At low 
temperature, the Luttinger liquids can be approximately 
thought of as an IEG consisting of quasiparticle excita- 
tions. Introducing mutual statistics or/ and Luttinger- 
type interactions among these excitations only shifts the 
value of X e f / . Thus the essence of Luttinger liquids is to 
have an IEG obeying FES as their fixed point. This is 
our characterization of Luttinger liquids in terms of FES. 

In this way, we have explicitly answered the three ques- 
tions raised in the introduction about Luttinger liquids: 

• The physical meaning of the Haldane's controlling 
parameter is the quasiparticle's effective statistics, 

X eff- 

• The Luttinger liquids, more precisely, the IEG, in- 
deed describe the infrared (or low-energy) fixed 
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points in 1-d systems, since their effective field the- 
ory at low energy is conformally invariant. How- 
ever, these fixed points are not isolated; they form 
a fixed-point line. Both the chemical potential and 
coupling constants are relevant perturbations that 
can drive the fixed point to move along the line, 
corresponding to the " renormalization" of the effec- 
tive statistics X e ff and leading to "non-universal" 
exponents. 

• It is conceivable that some strongly correlated sys- 
tems, exhibiting non-Fermi liquid behavior, in two 
or higher dimensions may also be characterized as 
having a GIG with appropriate statistics matrix 
as their low-energy or low-temperature fixed point. 
This is because the concept of exclusion statistics 
is independent of spatial dimeniosnality of the sys- 
tem. 

Moreover, we also showed that the effective field the- 
ory of 1-d IEG is a CFT with central charge c = 1 and 
compactified radius R = \/l/A. The particle-hole du- 
ality of the exclusons implies the CFT has an unusual 
duality R «-> 1/R, meaning that the CFT belongs to a 
new variant of the c = 1 CFT's, i.e the ones that are com- 
pactified on an "orbifold" S 1 fZi rather than on a circle. 
Physically, the differences are due to different constraints 
on the zero-mode quantum numbers. The CFT explana- 
tion makes a better understanding of the single-particle 
operators, especially, the anyonic (or exclusonic ) ones. 
Also, the CFT techniques provide a systematic way to 
calculate the correlation functions. 

Finally we observe several additional implications of 
this work: 1) Our bosonization and operator derivation 
of CFT at low energies or in low-T limit can be applied 
to Bethe ansatz solvable models, including the long-range 
(e.g., Calogero-Sutherland) one j|. 2) Here we have only 
consider one-species cases, i.e., with excitations having no 
internal quantum numbers such as spin. Our bosoniza- 
tion and characterization of Luttinger liquids are general- 
izable to GIG with multi-species, with the effective statis- 
tics matrix related to the dressed charge mat rix |9[ . 3) 
The chiral current algebra in eqs. (4.41) and (4.42) with 
A = 1/m coincides with that derived by Wen ]33[| for edge 
states in v — 1/m fractional quantum Hall fluids. So 
these edge states and their chiral Luttinger-liquid fixed 
points can be described in terms of chiral IEG. 

This work was supported in part by the U.S. NSF grant 
PHY-9309458, PHY-9970701 and NSF of China. 



right-left-moving representation. The density operator 
can be written as the Fourier transformations 

p(x) = p R (x) +p L {x), 

g>0 v 



q>a v 



; (e- iqx ~b q + e lqx bl), (Al) 



where Mr,l are given by M = Mr + Ml and b q = b- q 
for q > 0. 

The boson field <j>(x), which is conjugated to p{x) and 
satisfies 



[<f>{x),p{x')\ = iS(x-x'), 

is given by 

(x) = <Pr(x) + (/>l(x), 



(A2) 




(A3) 

with J = Jr + Jl- We have to assign the quantum num- 
bers such that there are only two independent variables 
in Mrl and Jr,l- A consistent choice is 



Then, 



Mr = J R , M L = -J L . 



J = Jr+Jl, m = j r -j l . 



(A4) 



(A5) 



Here </>o is an angular variable conjugated to M: 
[4>o,M] = i. The Hamiltonian ( p^ ) becomes 

H =\f dx[7rv N p(x) 2 + vj/7r (d x cj>(x)) 2 }, (A6) 
or by a field rescaling, 

H=^-£dx[U(x) 2 + (d x X(x))% (A7) 



where 



n(x) = ite-fpix), X(x) = e^^x). 



(A8) 



With X(x,t) = e iHt X(x)e~ iHt , the Lagrangian den- 
sity reads 



APPENDIX A: THE HARMONIC FLUID 
DESCRIPTION 

In coordinate space, there is a harmonic fluid descrip- 
tion |l(J of the Luttinger liquid. Instead of the 6-<j> rep- 
resentation that Haldane originally used, we prefer the 



C = -^d a X(x,t) d a X(x,t), 
2n 



(A9) 



which describes a free scalar field theory in 1 + 1- 
dimensions. 
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